Dissipated work and fluctuation relations for non-equilibrium single-electron 

transitions 
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We discuss a simple but experimentally realistic model system, a single-electron box (SEB), where 
common fluctuation relations can be tested for driven electronic transitions. We show analytically 
that when the electron system on the SEB island is overheated by the control parameter (gate 
voltage) drive, the common fluctuation relation (Jarzynski equality) is only approximately valid due 
to dissipated heat even when the system starts at thermal equilibrium and returns to it after the drive 
has been stopped. However, an integral fluctuation relation based on total entropy production works 
also in this situation. We perform extensive Monte Carlo simulations of single-electron transitions 
in the SEB setup and find good agreement with the theoretical predictions. 
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Statistical mechanics of small systems has been in the 
focus of intense interest over the past years. The common 
fluctuation relations, formulated, e.g., in Refs. PHI], m_ 
troduce equalities to describe irreversible processes. In 
the thermodynamic limit these are replaced by inequal- 
ities, of which the second law of thermodynamics is the 
best known one. The said equalities govern statistical av- 
erages over many repeated realizations of a given process 
driven by external control parameters. An individual re- 
alization is either dissipative or it extracts energy from 
the heat bath; however, the Jarzynski equality (JE) 



(1) 



should be valid Here (•} refers to averaging over an 
infinite number of repetitions, W is the work done in the 
driven process, AF is the free-energy difference between 
the equilibrium states of the system at the end points 
of the drive trajectory, and j3 = (ksT)^ 1 is the inverse 
temperature of the heat bath. The conditions for JE to 
be valid are amazingly few 0, foremost, the system 
needs to be in thermal equilibrium with the bath in the 
beginning of the drive. Typically, the validity and con- 
ditions of Eq. (JTJ have been theoretically discussed for 
abstract model systems [5, 6] , and experimentally for sys- 
tems, where controlled large sample averages are difficult 
to obtain [7|. 

A physical system of driven single-electron tunneling 
at low temperatures @ satisfies the requirements of con- 
crete experimental feasibility: averaging over large num- 
ber of realizations, and simple but accurate expressions of 
transition rates and energy relaxation become available. 
It has been recently shown theoretically [Io| and experi- 
mentally [ll[ that Eq. (JTJ) is valid for driven isothermal 
transitions in a single-electron box provided detailed bal- 
ance is obeyed. The general assertion is that this equal- 
ity is valid even if the system is overheated by the con- 
trol drive or it is driven to a full non-equilibrium state 
Here we demonstrate that this is not the case for 
the system we consider. Instead, we find that an inte- 



gral fluctuation theorem due to Seifert Q applies even 
in this situation. We expect the same conclusion to hold 
in general for driven, overheated systems where detailed 
balance is not obeyed. 

The single-electron box (SEB) jUGJ] (cf. Fig. [J) con- 
sidered in Refs. [1CH12| |. is a simple, yet a representa- 
tive system for our arguments. In a SEB, a tunnel con- 
tact admits electrons to enter or leave the island of the 
box. The electrostatic energy of the box with n (inte- 
ger number) excess electrons on the island is given by 
U(n,n g ) = E c {n - n g ) 2 . Here E c — e 2 /2Cs is the el- 
ementary charging energy of the box determined by the 
total capacitance Cs = C + C g + Co, where C is the 
capacitance of the tunnel junction, C g the gate capac- 
itance and Co the self-capacitance of the island. The 
relative energies of the different charge states n are de- 
termined by the gate voltage V g via n g = C g V g /e. In 
this work we discuss dynamics in the range < n g < 1, 
and at low temperatures, f3Ec 3> 1 such that the charge 
number can have values n = or n = 1 only. When an 
electron tunnels into (+) or out of (— ) the island, energy 
±[U(0,n g ) — [7(1, rig)] — ±2E c Sn g is released. We have 
written 5n g = n g — 1/2. The heat bath at temperature 
T is that of the phonons, and the lead of the SEB (on 
the left in Fig. [T] (a)) is assumed to be a reservoir in 
equilibrium with the phonon bath. However, unlike in 
Refs. [l(J [HI, nere we allow the small island of the SEB 
to be driven into non-equilibrium as a consequence of a 
gate protocol with many dissipative transitions (Fig. Q] 
(c)) and slow energy relaxation of the electrons. 

In Ref. [l(| it was concluded that Eq. ([T]) is always valid 
in a driven SEB if all the electrodes remain at the bath 
temperature T. This condition prevails if the drive in- 
jects non-equilibrium electrons (holes in out-tunneling) 
on the SEB island at a rate which is slower than the 
energy relaxation rate 7 of non-equilibrium excitations. 
The opposite limit / 3> 7, where / is the frequency of 
the cyclic gate drive, is the domain of the present dis- 
cussion. In this case there are two main regimes to con- 
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FIG. 1: Single-electron box (SEB). (a) Circuit diagram show- 
ing the electronic configuration and schematically the energy 
relaxation by wavy lines. Electrons tunnel at rates r be- 
tween the isothermal lead on the left and the small island in 
the center. The potential of the island is controlled capac- 
itively (capacitance C g ) by a variable gate voltage V g . The 
island is not necessarily isothermal due to the dissipative gate 
driving and weak energy relaxation (7) to the phonon bath. 
For more details, see text, (b) The energies of the two low- 
est lying charge states n = and n = 1 in the gate voltage 
range around Sn g = 0. The vertical red arrow depicts energy 
release in a n : — > 1 transition at a positive value of 5n g . (c) 
A section of the harmonic drive and a schematic presentation 
of the corresponding transitions of charge number n on the 
island, (d) Schematic presentation of the approximate island 
temperature evolution during the drive in quasi-equilibrium 
(regime (i), see text). 



sider as regards the energy distribution on the island, 
(i) Electron-electron relaxation rate is much faster than 
the drive, 7 C - C S> / S> 7, in which case the electrons 
on the island occupy a Fermi-Dirac distribution with 
higher temperature than the bath, see Fig. Q] (d). We 
emphasize that here we do not introduce two separate 
heat baths in the problem, but the second temperature 
is solely a result of the energy deposition on the island 
due to the drive by the control gate. If the protocol 
of the control parameter is stopped during the process 
to an arbitrary value, the system will relax to canoni- 
cal equilibrium corresponding to the value of the con- 
trol parameter. Thus the system obeys balance condi- 
tion in the sense of, and as required by the fluctuation 
theorems 0, [lj|. (ii) The injection rate is the largest 
frequency in the problem, / 3> 7c-c,7, which leads to 
a full non-equilibrium on the island, i.e., the electronic 
energy distribution deviates from the Fermi-Dirac one. 
If we consider the standard metallic SEBs, regime (i) is 
typical in experiments since 7 C _ C ~ 10 9 s _1 in an or- 
dinary metal, whereas 7 ~ 10 6 s" 1 for the correspond- 
ing electron-phonon relaxation at the sub-kelvin temper- 
atures where the SEBs are typically operated [TT1. [T6L Il7j| . 



In this work we discuss quantitatively the case (i) only. 
We consider a symmetric gate drive around the degener- 
acy, moving between charge states n = and n = 1 fl8j 
as depicted in Fig. [T](c). We assume unit peak-to-peak 
amplitude, where Sn g (t) varies between —1/2 and +1/2. 
We assume as usual that the system is in equilibrium 
with the bath before the gate drive starts at t = 0, as re- 
quested by the fluctuation relations in general. The drive 
ends after N back and forth ramps of the gate. The ther- 
modynamic work done is straightforward to obtain from 
its definition Q, or from the Markovian viewpoint jij. 
Since the beginning and end points of the drive are the 
same, Sn g = —1/2, there is no ambiguity of the proper 
expression of work to use [l9| , and the dissipation can be 
written as W — AF — J2i ±AUi, where the + sign refers 
to the rising n g half-periods, and — sign to the decending 
ones [Tol |. Moreover, for the sake of a simple argument, 
we assume in our analytical treatment that exactly one 
transition between the two charge states occurs in each 
half-cycle of the gate drive. This is a regime that can 
be achieved to a high accuracy - the probability of half- 
cycles with multiple transitions in a fully normal system 
to be considered below is « 3.0T 2 /(/3 4 £"^e 4 i?y), where 
T is the duration of the half-cycle and Rt is the junction 
resistance. Numerical results will be provided later for 
more general cases, where also multiple jumps in each 
leg can occur. For the ith half-cycle of the gate drive 
(increasing Sn g ), we may write in the single-jump ap- 
proximation by standard path-averaging arguments [l0( 



{e -PAU ) = 

T dTe-? Au ^e- % ^' r+ (-') r + ( T ) e - C 




(2) 



where AU = 2EcSn g is the work dissipated in the half- 
cycle, and Sn g is the gate position where the transition 
occurs. r + is the transition rate into the island that 
depends explicitly on time due to the gate drive, but 
also on the island temperature that can be different from 
that of the bath due to the dissipative transitions in the 
earlier half-cycles. The tunneling rates we apply in what 
follows are those for static biasing conditions, since we 
envision drive frequencies that are far below any relevant 
energy scale in the system. The descending half-cycles 
where an electron tunnels out of the island are identical 
in terms of the energetics and our argument. In what 
follows we assign the dissipated work in a multi-leg ramp 
as W — AF, and in a single leg as AU. 

If the island temperature stays constant at T, the de- 
tailed balance for the tunneling rates into and out of the 
island holds, r~(i) = e -^ AU ^T+(t), and by simple ar- 
guments presented in the supplementary on-line material 
we obtain from Eq. (fTT|) 
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Here Pr is the probability of making an: 1-^0 tran- 
sition in the corresponding reverse path, and the last 
step follows from our assumption of the symmetry of the 
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path and of exactly one transition in each leg. Since in 
the isothermal system, all the legs of the drive are inde- 
pendent, we see that JE is trivially satisfied in this case. 
This argument is more general in the isothermal process, 
not limited to the single-jump trajectories [Toj | . 

The situation is qualitatively different if the island is 
driven out of equilibrium. In what follows, we consider 
tunneling in a fully normal metal box, where 

r± = ^4fr / dE f( E ^)i 1 -fME±^U)} (4) 

are the tunneling rates when the lead is in equilibrium 
with Fermi-Dirac distribution f(E, /3) and island has 
a distribution /ne(-E)- For the full equilibrium case, 
/ne(£) = f(E,(3), we have T± = ±(e 2 R T )- 1 AU(l - 
e =F/3AC/ ) _1 . For the island in regime (i), /ne(-E) = 
f(E, /3 + A/3). A further simplified analysis in this regime 
can be carried out if we assume that the system stays 
close to equilibrium. We write A/3 for the instantaneous 
deviation of the inverse temperature of the island from 
that of the bath at the moment of the tunneling event, 
and assume that |A/3//3| <C 1. Linear expansion in A/3 
then yields 

r^r±(i-^r^), (5) 

We see immediately that detailed balance at temperature 
T is not obeyed any more if A/3 ^ 0, but 

T-/T+ ~ e-P AU [l - -At/A/3]. (6) 

We find a corresponding linear correction to the expres- 
sion in Eq. (| 1 1 1) for the ith half period of the gate drive 
as 

(e^ AU )cl-l(Q)Af3. (7) 

Here (Q) is the average value of heat dissipated in cross- 
ing the degeneracy, which is evaluated for the isothermal 
system here: according to Ref. [l(|, (Q) > is approxi- 
mately proportional to the sweep rate. 

Although (e~ /3AC/ ) ^ 1 in Eq. in general, this of 
course does not prove JE wrong, since the requirement 
of equilibrium at the beginning of the gate (half-)period 
was not imposed here. However, we may make use of 
Eq. ([7]) to prove that JE does not hold in the following 
simple example where the ramp starts under the equilib- 
rium conditions. Suppose that the symmetric gate ramp 
consists of only two linear legs (N = 1), see inset in 
Fig. [3] Again exactly one tunneling event occurs in each 
period. Next, we assume that the first tunneling event 
heats or cools the island by an amount determined by its 
heat capacity C, such that AT = before the first tun- 
neling, and AT — E/C after the first tunneling, where 
AT = — fc^ 1 A/3//3 2 and E is the energy deposited on 
the island by the tunneling electron. Furthermore, we 
assume that the relaxation of heat is so slow that the 



temperature of the island is changing only at the transi- 
tions during the sweep (20j . In this sweep, we may then 
write similar evolution for the full cycle as in Eq. (jlll) . 
but now assuming that the temperature on the island de- 
pends on the energy at which the electron tunnels in the 
first leg. This analysis is presented in the supplementary 
on-line material. After straightforward algebra we obtain 
for the dissipated work W - AF = AUi - AU 2 

(e-^-^))^l-^/3 2 (Q) 2 . (8) 

This equation shows that although (e~^ AUl ) = 1 for 
the first half-period with A/3 = 0, and (e" /3AC/2 ) > 1 for 
the second half-period with typically A/3 < 0, the over- 
all average assumes values ( e -0( w - AF )} < i j d U e to the 
heating induced correlation between the two legs. We 
note by using Jensen's inequality, < (e~ x ), that 

this conclusion is consistent with the second law of ther- 
modynamics. 

To check our analytic predictions we performed numer- 
ical simulations of single electron transitions in the SEB 
setup by using the standard stochastic Monte Carlo sim- 
ulation method. Unlike in the analysis above, in these 
simulations no approximations were made: the number 
of jumps within a leg was not restricted to one, and we 
used the exact form of the tunneling rates instead of the 
linearized approximation of Eq. (1151) . Furthermore, in 
the analytical results, we have expanded the tunneling 
rates up to the second order in A/3 for more precise 
comparison. We set the resolution of AU(t) G [—1,1] 
to 0.0005 using linear increments and the resolution of 
the temperature difference between the lead and the is- 
land AT e [0,2T] to 0.001T. The control parameter 
protocol consisted of a linearly increasing half sweep 
n g (t) : [0,tf],n g (t) — t/tf and a linearly decreasing one 
n g (t) : [tf,2t f ],n g (t) = 2-t/(2t f ), with 4000 steps. We 
set the sweep rate / = 2/tj to 10 7 s _1 , which gives about 
95% single jump legs. In each simulation the tempera- 
ture of the heat bath was set to T = 0.1 K, the charging 
energy to Eq = 40fceT and the tunneling resistance to 
R T = 450 kfl. 

We calculated the tunneling rates (Eq. (fT3"|) ) numer- 
ically using the standard trapezoidal rule integration. 
We determined the integration limits and the number of 
trapezoids by doubling both of them until the error was 
less than 10~ 15 for all values of AU and A/3. The heat 
generation method is presented in detail in the supple- 
mentary material. We assumed weak coupling between 
the electrons on the island and the phonon bath. Thus 
the temperature of the island was controlled only by the 
heat generation due to the tunneling events during the 
sweep. 

To numerically test the prediction of Eq. (8) for 
^g-^(VK-AF)^ we p er f ormec i i 75 x 10 11 repetitions for 

different values of the heat capacity C ranging from 
ll-Ec/K to SO-Ec/K- The simulation data are shown 
in Fig. [2J First, as expected JE is accurately satisfied in 
the limit of AT = 0. This was verified with several differ- 
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FIG. 2: Simulation results (triangles) of (e" fi(w '- Af) ) and the 
corresponding theoretical approximation of the single jump 
and small AT/T approximation, Eq. (8) (dashed line). The 
horizontal axis shows the average temperature after the first 
half sweep. The corresponding values of the heat capacities 
from left to right are: 50E C /K, 30E C /K, 20E C /K, 15E C /K, 
12Ec/K and llEc /K. The error bars are the standard er- 
ror of the mean of the corresponding data. Each data point 
is obtained from 1.75 x 10 11 independent realizations of the 
tunneling process. The insets (a) and (b) show the sampling 
from simulation with C = 12E C /K with all the 1.75 x 10 11 
repetitions, (b) demonstrating explicitly that even the tails of 
the distribution are well sampled. 



ent sweep rates (data not shown). Most importantly, we 
find that (e~^^ w ~ AF ^) decreases with increased heating 
in good agreement with the theoretical prediction [21|. 

To explicitly demonstrate the details behind the vio- 
lation of the JE, we calculated (e _/3AC/ ) separately for 
the increasing and decreasing parts of the sweep. The 
results are shown in Fig. [3l The average is very close 
to unity for the first half of the sweep, which starts in 
equilibrium. The slow increase is a result of the contri- 
bution from multijump legs. However, during the second 
half sweep the initial temperature of the island is higher 
than that of the bath and we obtain (e-P AU ) > 1 for 
that leg. This result is in agreement with the theoretical 
prediction of Eq. (7). 

An integrated fluctuation theorem (IFT) based on to- 
tal entropy production, including the entropy change of 
the system and the surrounding medium along a single 
trajectory, Astot, was discussed by Seifert in [8j]. It was 
concluded that IFT in the form 



-Astot 



1 



(9) 



applies without either the assumption of equilibrium 
starting condition or detailed balance. We can apply 
the single-trajectory entropy production, with the as- 
sumptions in the analytic treatment of this paper in a 
straightforward manner, and obtain the relation 
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FIG. 3: Simulation results of {e~ f) ^ u ) over the first and the 
second half-sweep as a function of the scaled average tem- 
perature after the first half sweep. The corresponding values 
of the heat capacities for the data points are the same as in 
Fig. [2] The starting point of the second sweep is the tem- 
perature difference produced in the corresponding first half 
sweep, while the first half sweep starts from equilibrium. The 
statistical errors are of the size of the data points. Each data 
point is obtained from 1.75 x 10 11 independent realizations 
of the tunneling process. The inset (a) illustrates the control 
parameter protocol during the first half and the second half 
of the sweep. 



By applying the single-jump techniques above, we find 
that (A[/e~ ,9AC/ ) = —(Q), and inserting this result and 
Eq. O into Eq. (flu), we conclude that Eq. © is valid 
identically (at least within the approximations made) 
for a leg with arbitrary island temperature. Therefore 
Eq. ([9]) is valid also for the overheated two-leg trajectory 
where Eq. JT]) fails. 

In summary, we have analyzed the prediction of the 
Jarzynski equality for driven single-electron transitions 
under experimentally relevant conditions: the system 
may be either in equilibrium, in quasi-equilibrium (over- 
heating) or in true non-equilibrium. The JE is satisfied in 
the first case. We have shown both analytically and nu- 
merically that Eq. (JTJ) is not applicable when dissipative 
heat causes a temperature change in the SEB, although 
the system starts at equilibrium and returns to it. How- 
ever, the integrated fluctuation theorem of Eq. ([9]) holds 
even in this case. The experimental realization of our 
setup requires low tunneling rates and slow relaxation 
rate between the system and the bath: these require- 
ments are best satisfied using tunnel junctions between 
normal metals and superconductors [17j . or in semicon- 
ducting quantum dots [12, HH . 
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I 

I. DERIVATION OF EQ. (3) OF THE MAIN TEXT 

We start with Eq. (2) of the main text: 

( e -/3A£/ )= f T dTe -^m-r) e -£dT>r+(T') r +( T)e -jTdT'r-(T>)_ (u) 
Jo 

Inserting the detailed balance condition T~(t) — e~P Au ^T + (t) into this equation we may rewrite it into 

{e -PAU )= f T ^ dre -/o r ^T+(r') r - (r)e -^dr'r-(r'). (12) 
JO 

Next we make a change of variable in the main integral, t = T — r, and note that = T^(T — u) for the reverse 

trajectory R and for any time instant u along the symmetric trajectories considered here, yielding 

(e-^ u ) = [ T dte-tf-^-niT-T'^-^-f^dT'r-iT-T')^ (13) 
Jo 

Making the change of integration variable in the integrals of the exponents as well, r" = T — r', and reordering the 
terms in the integrand, we finally obtain 

( e -f>*V)= [ T dte-ti dT " r n^r-(t)e-S? dT '' r x^''\ (14) 
Jo 

By direct inspection we identify this as the probability of making exactly one n : 1 — > transition in the reverse 
trajectory, i.e., Pr, and due to the condition that exactly one transition is taking place in its mirror trajectory, we 
finally conclude that 

{e -PAu ) =Pr = 1 (15) 

in this case. 



II. DERIVATION OF EQ. (7) OF THE MAIN TEXT 

We first expand the tunneling rate into the box 

r+ = / dE K E > m-f{E + AU, 8 + Aft)] (16) 

for a small inverse temperature difference A3 up to the linear correction. Since df(E, (3)/d(3 = —Ef(E, /3)[1 — f(E, j3)], 
we obtain 

^ + -^ J dEf{E,l3)[l- f{E + AU,f3)] + -^L J dE{E + AU)f(E, B)f(E + AU,B)[\ - f(E + AU,B)}. (17) 

Since all the temperature arguments in Fermi distributions are now equal (= B), we may drop them for now. Equation 
(jTTJ) can thus be written as 

T+ ~ J dEf(E)[l - f(E + AU)} + J dE(E + AU)f(E)f(E + AU)[1 - f(E + AU)}. (18) 

Since f(E)[l - f(E + x)] = [f{E] - f(E + x)}/(l - e~^ x ), we obtain from Eq. (fLS]) 



where 



r+ ~ T+ + ^ 1 _ ^ pAU J dE(E + AU){f(E) f(E + AU)}f(E + AU), (19) 

r ° = dEf{E,B)[l-f{E + AU,B)] = i (20) 

is the tunneling rate into the box at equilibrium temperature. Let us derive the analytic expression of the integral in 
Eq. ([T9]) step by step for illustration. 

J dE(E + AU){f(E) - f(E + AU)]f(E + AU) = 

J dE(E + AU)f(E + AU)[1 - f(E + AU)} - J dE(E + AU)f(E + AU)[1 - f(E)} = 

J dE' E'f(E')[l - f(E')] - J dE' E'f(E')[l - f(E' - AU)}. (21) 

Using 

J dEEf(E) [1 - f(E + x)} = -~ 1 _ X e _ px , (22) 

we obtain from Eq. (|2"Tj) . 

J dE{E + AU)[f(E) - f(E + AU)]f(E + AU) = ± - A ^ AU . (23) 

Combining Eqs. (fT9|) and ff23|) . we get 

, , AB AU AU e 2 R T 

where 

r " = ^ JdEf {m - f(E - AU)] = (25) 

is the tunneling rate out from the box at equilibrium temperature. 
Similarly we find for the opposite rate 

r-^r -(i-^r+A/?). (26) 
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Therefore, instead of the original detailed balance Tq /T = e^ AC/ , we find 

^ ~ + ^(r+ - r )A/3] = e^ u [l + IaUAP], (27) 

which is the reason why JE fails. Here, in the last step we have used Tq — Tq = A[//(e 2 i?T)- 

Now consider a symmetric trajectory of the gate from time to r. We assume that one and only one jump occurs 
when crossing the energy degeneracy. We find for this leg 



( e -/3AE/ )= / e ^AC/(r) e -/; rf r'r+(r') r + (T)e -/ T r rfrT-(r') rfr 



T e-f Td -T + (T')r-(T)e-{?dT'r- ( r') dT+ l A(3 C Af/(r)e - / T ^T + (r') r - (r)e - /* dr'T-(r') dT _ (2g) 

2 .In 



In the second step we have made use of the result Eq. (|27| . Note that from now on we can use the expressions or 
values of the quantities at the initial temperature, since the corrections would yield errors in the next order only. 
Now, we find that the second line of Eq. (|2"5)l can be rewritten such that 

(e-P AU )~P R +~(Q) R Al3, (29) 

where P R is the probability of the transition in the reverse trajectory over the same gate section, and (AU)r is the 
expectation value of AU over this reversed trajectory. Since we assume linear trajectory and the assumption of success 
with one jump each time P = 1, we have P R = P =1 and (Q)r = — (Q). Thus, 

(e-^ u )^\- l -(Q)Ap. (30) 



We can write (l30l also in the form 



(e-^)^l-i<Q)A/3. (31) 



Since (Q) > 0, and A/3 < 0, we find that (e l3AU ) > 1 in the case of an overheated island. 



III. DERIVATION OF EQ. (8) OF THE MAIN TEXT 

Now we consider a back-and-forth gate trajectory as in the main text. We make use of Eq. (f3"T|) 

(e-? AU )^l~ l -(Q) A/3 (32) 

to prove that JE is not satisfied. Here j3 = l/fc^T, where T is the bath temperature. We consider a simple trajectory, 
which consists of a linear n g : — > 1 sweep followed immediately by a similar but opposite sweep n g : 1 — > 0. We 
make the simplifying assumption that in each half-period, one electron tunnels in the preferred direction. We also 
assume that when heat is released in the island, it adjusts its temperature instantaneously, but the relaxation rate of 
heat to environment is much longer than the duration of the sweep. In this sweep, lasting over a period 2r = r + r, 
we may then write the average as 

/oo pIt pT 

dE dr 2 drie-' 3 [ AC/ ^)- At/ ^)] e -io T1 r + (A£/(r'),To)dr' 7+(E)ri)x 
-oo Jr JO 

e -f^r-(AU(r'),nE))dr' T _r AU{T2)>T{E) y (33) 

Unlike in the previous se section, we have dropped in the integrals the probability of the second transition, since we 
explicitly assume that only one transition occurs in each leg. Here, T±(AU(t),T') are the tunneling rates into/out of 
the island when the island has temperature T", and 

1+ {E, Tl ) = -^-f(E - AU{ Tl ),m - f(E,j3)] (34) 
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is the corresponding equilibrium transition rate density of r + at energy E. The expression in Eq. ([33| assumes that 
the system starts in equilibrium (temperature T), and the tunneling event in the first leg influences the temperature 
via T(E) = T + AT(E), and therefore the tunneling probability out from the box in the second half of the cycle. We 
assume that 

AT(E) = |, (35) 

i.e. the increase of temperature is proportional to the energy E deposited on the island, and C is the heat capacity. 
Since 

T_(AU,T') = J dEf(E,l3')[l-f(E-AU,0)] = -^- J dEf(E + AU, /3)[1 - f(E, /?')] = T + (—AU, T')(36) 



always, we may write formally Eq. f)33[) as averaging over two n g : — > 1 sweeps at the two different temperatures as 

/OO pT 
dE / dne-^^e-^r+^^^'j+iE,^) x 
- oo JO 

dT2e -/3A C /(r 2 ) e -^r + (A[/(r'),T( £ ))^ r+ (A[/( T2 ),T(i?)). (37) 



Here the second line is the average of Eq. ([32]) with A/3 = -AT(E)/(k B T 2 ). Thus we may write (J3TJ) as 

dE r dTie ^ A ^)e-^ 1 r + (A^^T) d r' 7+( ^ Ti)(1 + ^^M i?) 
Jo 2 C 

drie -PAU(T l)e -!? r + (Aa(r'),T) d r' r+(At/(ri))T) 

2 C j 

[ |°° dE7£7/(£-A^(ri),T)[l-/(B,T)]] =1+^^M<AI7) K . (38) 
Here, (AU)r — —(Q) is the average value of AU in the reverse trajectory, and finally 

{e -P(W-AF) ) ^ 1 _l^p2 {Q) 2_ m 



